ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ ЗНАЧЕНИЙ ВРЕМЕННЫХ РЯДОВ 
МЕТОДОМ ДИНАМИЧЕСКИХ УСЛОВНЫХ КОРРЕЛЯЦИЙ НА ОСНОВЕ 
МНОГОМЕРНОГО АСИММЕТРИЧНОГО РАСПРЕДЕЛЕНИЯ ЛАПЛАСА 

О.А. Бельснер, О.Л. Крицкий 

Томский политехнический университет 
Е-таіІ: оІедкоІ@трН.рНШри.есіи.ги, ЬеІ5пег@Ьк.ги 

Предлагается модификация модели динамических условных корреляций. Модификация заключалась в отказе от предположе- 
ния о многомерном нормальном распределении логарифмов приращений дневных котировок финансовых инструментов, 
вследствие ненулевых значений коэффициентов асимметрии и эксцесса, и в использовании асимметричного многомерного ра- 
спределения Лапласа, позволяющего моделировать линейные комбинации одномерных случайных величин, что особенно важ- 
но при расчете предельной величины риска ( Ѵаіие-аі-гізк , ѴаК) для портфелей финансовых инструментов. 


Понимание механизма взаимодействия финан- 
совых величин является проблемой, довольно дли- 
тельный период рассматриваемой в экономике. 
Переменные оказывают влияние друг на друга не 
только через первый, но и через второй моменты 
совместных распределений. Это означает, что из- 
менение значения одной переменной может воз- 
действовать не только на уровень другой, но и на 
степень изменчивости остальных переменных. Для 
достижения состоятельной, несмещенной и эф- 
фективной оценки данной зависимости предлага- 
ется большое количество эконометрических мето- 
дов, позволяющих осуществить наиболее точный 
анализ. 

В течение последних двадцати лет в финансо- 
вой эконометрике происходит активное развитие 
моделей, описывающих процессы изменения цен, 
основной особенностью которых является непо- 
стоянство безусловной или условной дисперсии. 

Причины изменения дисперсии у финансовых 
показателей различны. Среди прочих можно выде- 
лить непредвиденные политические события, рез- 


кие колебания уровня предложения на рынке (на- 
пример, вследствие поступления на рынок больших 
объемов валюты, товаров или ценных бумаг) или 
его сокращение (что, согласно закону спроса и 
предложения, приводит к росту цен). Очевидно, что 
изменения дисперсии цен финансовых инструмен- 
тов оказывают существенное влияние на финансо- 
вые сделки вследствие увеличения рисков потерь. 

Общий подход к построению моделей с изме- 
няющейся дисперсией предполагает, что значение 
финансового показателя X, в момент времени ( 
определяется посредством уравнения: 

X, = ц(() + о,е І , 

где ц(і) - условное математическое ожидание про- 
цесса Х„ ь',~УѴ(0,1) - стандартная нормально ра- 
спределенная случайная величина, сг, - условное 
стандартное отклонение. 

Изучение изменений волатильности привело к 
появлению в начале 1980-х гг. класса моделей авто- 
регрессии условной гетероскедастичности (АКСН) 
[1] вида: 
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а; =сс + ^Р : е 2 ІЧ , 

І = 1 

где а, Д, і=1,р - некоторые скалярные параметры, 
причем а>0, Д>0. В АКСН-моделях предполагает- 
ся зависимость волатильности от множественных 
шумов вида е, 2 ~А(а„сг, 2 ). 

За последние двадцать лет были разработаны 
многочисленные модификации АЯСН-моделей, 
позволяющие отказаться от предположений о неза- 
висимости волатильности от своих предыдущих 
значений и учесть автокорреляцию в них. Наибо- 
лее известной модификацией АКСН-модели явля- 
ется обобщенная модель авторегрессии условной 
гетероскедастичности (САКСН) [2] вида: 

= а + І + Ъі*м ’ 

м І = і 

где а, Д, у р і=\,р,]=\,д коэффициенты а>0, Д>0, 

ъ> о. «+ХД+Ха < 1 - 

і і 

САКСН позволяет моделировать изменяющую- 
ся во времени условную дисперсию в виде регрес- 
сионной зависимости от предыдущих ее значений 
по времени. Модели типа САКСН находят широ- 
кое применение в различных сферах, однако наи- 
большее распространение они получили в финан- 
совой эконометрике при исследовании финансо- 
вых рынков, так как их существенным преимуще- 
ством являются свойства быстрого реагирования 
на любые наблюдаемые изменения и восстановле- 
ния после сильных колебаний рынка [2]. 

В настоящее время имеется большое количе- 
ство модификаций методов АКСН и САКСН. Наи- 
больший интерес представляют многомерные 
обобщенные модели авторегрессии условной гете- 
роскедастичности (МО АКСН), учитывающие за- 
висимость волатильности. Задачи оценки стоимо- 
сти активов, опционов, выбора портфеля и расчета 
величины ѴаК являются примерами успешного 
применения моделей М САКСН [3]. 

Модели МО АКСН могут быть использованы 
для оценки степени влияния волатильности фи- 
нансовых рынков на такие величины как экспорт и 
производительность, а также при расчете изменяю- 
щихся во времени хеджирующих соотношений. 

Основным недостатком использования МОАКСН 
является то, что число оцениваемых в них параметров 
резко увеличивается по мере роста числа переменных. 
В большинстве случаев уже при 4-х рассматриваемых 
переменных произвести оценку соответствующих па- 
раметров затруднительно (минимальное число для 
различных моделей МОАКСН равно 20). 

В 1990 г. Т. Воііегеіеѵ [4] предложил класс моделей 
с постоянными условными корреляциями, предпо- 
лагающий, что условные ковариации пропорцио- 
нальны произведениям соответствующих условных 
стандартных отклонений. Данное ограничение зна- 
чительно уменьшает число неизвестных параметров 


и упрощает оценку. Тем не менее, предположение о 
постоянстве условных корреляций на практике ока- 
зывается недостижимым. В данной работе предлага- 
ется к рассмотрению класс моделей динамических 
условных корреляций (Бупатіс Сопбіііопаі Соггеіа- 
1іоп8, БСС [5]), учитывающий нелинейную природу 
изменения финансовых временных рядов и дина- 
мические свойства корреляций. 

Модель динамических условных корреляций 
была предложена К. Егшіе и К. ЗНеррагсІ в 2002 г. [5]. 

Рассмотрим некоторый «-мерный процесс для 
доходности некоторого актива Х,еК т ' п , /=Т...Г[5]: 

х, = Иі (Ѳ)+нГ(Ѳ)- 8 ,, 

И,(Ѳ) = Е(Х, |Ц_,), 

И, = Ѵаг(Х ( |Ц_,), (1) 

где О, - фильтр (доступная информация) на мо- 
мент времени г, а е, — нормально распределенные 
величины. 

Пусть Н, - функция одномерных динамических 
переменных и динамических линейных корреля- 
ций вида: 

н,=ОАЦ> (2) 


Тогда 0, можно определить по методологии 
САКСН(Ь,8): 


а= д 


1=1 5=1 

й* = (1іа ё||Ѵ^'||. ( 3 > 

где Д - диагональная матрица с элементами лД 
і=\...п, 1=1. ..Т, 0 - безусловная ковариационная 
матрица е„ 8=0, 'г.. Предполагается, что е~ Щ ОД), 
0=Е(е,е,), а„ Д - параметры (скаляр), удовлетво- 


і X 

ряющие условию: ^а, + ^ Д <1, К, - корреля- 

1=1 5=1 

ционная матрица, Д - положительно определен- 
ная матрица, характеризующая волатильность. 

Запись К, в виде (2) гарантирует универсаль- 
ность корреляционной матрицы, а выражения для 
записи Д и 0, обеспечивают положительную опре- 
деленность д. 

Оценка параметров в (3) осуществляется двух- 
шаговым методом максимального правдоподобия. 
При этом логарифмическая функция правдоподо- 
бия разбивается на две: одна для параметров, опи- 
сывающих одномерную волатильность, другая - 
для параметров, определяющих корреляцию (так 
называемый двухшаговый БСС метод). Такой под- 
ход позволяет оценить системы, содержащие боль- 
шое количество параметров. 
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На первом шаге вектор Ѳ разбивается на два 
подвектора Ѳ=(^,<р), где ^ содержит параметры во- 
латильности, а (р — корреляции. Тогда логарифми- 
ческая функция правдоподобия запишется в виде: 

Ь т (Ѳ) = І т ^) + Ь г {ф. 

Следовательно, процедура оценки на первом 
этапе сводится к определению максимума функ- 
ции На втором шаге, после нахождения 

оценки переходят к максимизации функции 
Ьт((р\^) и нахождению оценки ф, принимая во вни- 
мание, что г,=Д “V,. 

Для доходов финансовых активов свойственны 
ненулевые значения асимметрии и эксцесса, 
вследствие этого их нельзя описать с помощью 
нормального распределения. Выделяется несколь- 
ко видов многомерных распределений, характери- 
зующихся высоким уровнем коэффициентов эк- 
сцесса и асимметрии. Однако большинство таких 
распределений либо является слишком сложным 
для оценивания при помощи метода ОАК.СН, либо 
обладает нежелательными свойствами (например, 
бесконечная дисперсия), что, естественно, снижа- 
ет степень их использования для описания финан- 
совых процессов. 

Одним из распределений, достаточно точно 
описывающих основные свойства доходов финан- 
совых инструментов и обладающих определенной 
гибкостью, является асимметричное многомерное 
распределение Лапласа, преложенное 8. Коіг, 
Т..І. КогиЬошкі и К. РоскопАі в 2003 г. [6]. Одно- 
мерное распределение Лапласа или двойное экспо- 
ненциальное распределение широко используется 
в финансовом моделировании. В [6] асимметрич- 
ное многомерное распределение рассматривается 
как подкласс геометрически устойчивых распреде- 
лений (устойчивых при геометрическом суммиро- 
вании). Данное свойство позволяет использовать 
асимметричное многомерное распределение Ла- 
пласа для моделирования линейных комбинаций 
случайных величин, подчиняющихся одномерному 
симметричному распределения Лапласа, что явля- 
ется особенно важным при расчете коэффициента 
ѴаК для портфелей финансовых инструментов. От- 
мечается, что данная особенность ранее считалась 
характерной лишь для Парето-устойчивых распре- 
делений и их частного случая, используемого наи- 
более часто - нормального распределения. Харак- 
терная черта класса устойчивых распределений со- 
стоит в том, что сумма случайных величин, подчи- 
няющихся какому-либо устойчивому распределе- 
нию, также распределена по этому закону. Устой- 
чивые распределения, как правило, отличаются от 
нормального более медленным приближением 
плотностей распределений к оси ординат. Иными 
словами, на их «хвостах» накапливается вероят- 
ность более значительная, чем на «хвостах» нор- 
мальных распределений с теми же параметрами. 

В геометрически стабильной модели [6] доход 
г м рассматривается как сумма меньших доходов і м> 
в течение периода времени Лр). 


Лр) является случайной величиной с геометри- 
ческой функцией распределения вероятности 
Р(Лр) = Л = р( 1 —р) н ' , к= 1,2,.... Геометрически ста- 
бильное распределение может быть приближено к 
нормированной сумме геометрически стабильных 
моделей, если параметр р случайной величины Лр) 
стремиться к нулю. Иными словами, некий массив 
случайных величин ^подчинен геометрически ста- 
бильному распределению на множестве 9?" тогда и 
только тогда, когда [6]: 

Яр) 

а{р)'^ і {к{р) + г и) ) — при р >0, (4) 

(=1 

где {г м =(г 1 ( '' ) ,...г„ ( ‘ л ), д>\) - последовательность неза- 
висимых одинаково распределенных случайных 
векторов в пространстве 9?", не зависящих от Лр), 
а(р)> 0 , к(р)е)Я ", — означает сходимость по ра- 
спределению. 

В данном случае многомерное асимметричное 
распределение возникает, когда существует огра- 
ничение для предела распределения - конечность 
второго момента. 

В [6] было показано, что в случае, если каждый 
вектор г характеризуется средним т„ і=\,п и дис- 
персией сТу, /=1,и,у-1,и, а а(р)=Фр и к(р)=т(Фр-\), 
случайная величина X, определяемая выражением 
(4), подчиняется АМБ-распределению с характе- 
ристической функцией: 

4X0= , ' , 

1 + -('Н(-і('т 
2 


где /е9?" х " и Не\ К" х “ - положительно определенная 
матрица с элементами /=/ на главной диагонали 
и недиагональными элементами а р /У/. 

Функция многомерного асимметричного ра- 
спределения Лапласа (АМБ) для изменяющихся во 
времени Н , и г, имеет вид: 


ДО 


2ехрДЯ,Ли) Г г/дД Ч 
(2тг)" ,2 \Н,\ т [2 + т'н; 1 гп; 



где ѵ=(2-п)/2 и К ѵ (и) - модифицированная функ- 
ция Бесселя третьего рода, 


. _ ( И /2УТ(1/2) | , т 

Г(ѵ +1/2) ] 


и> 0, ѵ>-1/2. 


Здесь вектор т - параметр положения, а матри- 
ца Н - коэффициент масштаба данного распреде- 
ления. 

Одной из важных особенностей АМБ-распреде- 
ления является тот факт, что оно унимодально с ну- 
левой модой. Вследствие данного факта параметр 
т не только определяет среднее данного распреде- 
ления, но также характеризует уровень асимме- 
трии. 
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Обобщим предложенную модель динамических 
условных корреляций, для чего предположим, что 
остатки подчиняются АМЬ-распределению. 

Логарифмическая функция правдоподобия в 
случае асимметричного многомерного распределе- 
ния Лапласа с учетом того, что Я,=ДЛД, имеет вид: 

(Ѳ) = X (Г' О-'К'Ц-'т - \ ІП \Б , Я , Б , | + 

+^1п(г;д- 1 д-ДД )-ЩтЪ ; 1 к ;'о;'т) + 

+ НК ѵ ф + тЪ ; 1 К , 7) »(/;Д - 

— ~ 1п(2тг)). 

Применим двухшаговую процедуру оценки: 

1 аш . (Л) = X и |ду. | + 

+ 2 1п (л А Ул ) - М'и Д" 2 Л + 

+ Іп(/С ѵ ( л /(2 + т’0; 2 І к т)( Г ;о; 2 І к г,)) - " 1п(2тг)). 


Так как е,=Д V, и е',=т Д ', имеем: 


АмА^)=Х(^УО~ЦН*'1 + 

г=1 ^ 

+ ^ 1п(г ( ' д _1 д _1 дД ) - 1п(/и Д ; 'я о ; ] т) + 
+ЩК ѵ ф + ш'д-'/дД»(/-Д ,-Д Г Д ; V ( )) - 
+ ~ 1п(27Г )). 


Если принять й= 25+3, где 5=0, 1,2,..., то 

5 1 + 1 

ѵ = - — , К=К_ Ѵ . В данном случае функция Бес- 
селя может быть записана в виде: 


т= 

( г ( Н~ 1 т — 


уА + т п ( т) ехр| 

А 

^(2 + тН 1 х т){г 1 Н 1 х г 1 ) ^ 

(2тг у] 

'РА 

X 

н, 

1/2 


, (■? + *)! 

х Х {$-к)\к\ (2уІ(2 + т'Н;'т)(г;н;\)Т к . 


Следовательно, логарифмическая функция 
правдоподобия примет вид: 


= Х^1ч(2+т'о;7гД; 1 т) +г/д-Д-Д> - 

-р + т 1 д-‘< ДЖДДд-Д,-^ Ч) + 

іп{У ( - ; + А ~ )! X 
(з-к)\к\ 

х(2р + т 'д-Д,- ‘д-'шіігХ'д- Д; V,))-*} - 

-(5 + 1)1п(2л:) - 

-^ІпСд/г/Д-Ч-ДЛ ) ■ — |іп |ДД Д |}. 
Применим двухшаговую процедуру оценки: 
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Ам^) = Х^ 1п 2 + Х? 


V 


^14, 


2Л м/ 


.2 X 


X- 


т 2 У^ А > 




І=1 


V "Й 


2 ѴѴ 


ІП 


(5 + Аг)! 


х 

І? 0 (8-к)\к\ 


Ні у Ѵ'= 1 "И ) 
\~ к 


" т 2 Л( " г 2 7 

2 + Ѵ^- Ѵ^_ 
ыКк&Ьи; 


-С?+1)1п(2тг). 


-^іп( г ;д-Ч')-ііп|д І+^'Д-ЧО' - 
-р + е ;к;'(е;)')(еХ'е, ) + 
-^ІХ гЧѵД , (2л/(2 + <Д"«)ЖУЧ )Г*| - 

-С?+1)1п(2тг)}. 


Применяя метод максимального правдоподо- 
бия, получим систему: 

'аУмД|) = 0 

дК 

^АмЛ ( Р\Л) _ - 

ад ' 1 

Таким образом, предложена модификация мо- 
дели динамических условных корреляций (2). Так 
как предположение о нормальном распределении 
не отражает действительного поведения рисковых 
активов на финансовых рынках, было использова- 
но многомерное несиметричное распределение Ла- 
пласа. Получена система уравнений (5), решение 
которой однозначно определяет неизвестные коэф- 
фициенты в (3). Найденные коэффициенты детер- 
минируют предлагаемую модель и позволяют ис- 
пользовать ее при анализе эмпирических данных. 
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